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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


if 

i 

{ 


U.  S.  customary  units  of  measurement  used  in  this 
to  metric  (SI)  units  as  follows: 


_ Multiply _  _ By 

cubic  inches  16.38706 

inches  25.4 

inches  per  second  25.4 

miles  per  hour  (U.  S.  1.609344 

statute) 

square  inches  6.4516 

square  inches  per  6.4516 

second 


report  can  be  converted 

_ To  Obtain _ 

cubic  centimetres 
millimetres 
millimetres  per  second 
kilometres  per  hour 

square  centimetres 

square  centimetres  per 
second 
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FREQUENCY  RESPONSE  CHARACTERISTICS  OF  MILITARY  VEHICLES 

PART  I:  INTRODUCTION 
Background 

1.  Driver  safety  standards  and  the  design  specifications  of  new 

electronic  and  optical  equipment  for  vehicles  require  estimates  of  the 

values  of  the  acceleration  components  expected  to  occur  at  a  point  on  a 

vehicle  when  it  crosses  an  obstacle  or  other  terrain  feature  at  a  given 

speed.  The  electronic  and  optical  equipment  have  shock  and  vibration 

criteria  that  are  expressed  by  maximum  acceleration  tolerances.  The 

driver  tolerance  is  generally  expressed  as  a  limit  on  the  power  absorbed 

1  2 

which  in  turn  can  be  related  to  acceleration.  ’  These  criteria  limit 
the  speed  at  which  a  vehicle  can  cross  obstacles  and  other  terrain 
features. 

2.  A  dynamical  model  of  a  vehicle  is  generally  required  to  predict 

the  acceleration  components  at  a  point  on  a  vehicle  as  it  runs  over  a 

specified  terrain.  This  model  consists  of  a  combination  of  spring-mass- 

damper  elements  that  are  assembled  to  approximate  as  closely  as  possible 

the  measured  dynamic  response  of  wheel  and  track  laying  vehicles.  Many 

models  of  this  type  have  been  developed  and  used  with  varying  degrees  of 
3  4 

success.  *  Within  the  limits  of  this  formalism  the  predicted  dynamic 
response  of  a  vehicle  will  always  be  model  dependent.  In  order  to  tepre- 
sent  as  closely  as  possible  the  actual  interaction  of  a  vehicle  with 
terrain  features  it  is  of  value  to  develop  a  model  independent  methbd  for 
predicting  the  dynamical  response  of  vehicles.  Such  a  method  can  be 
developed  to  predict  the  root  mean  square  (RMS)  values  of  the  acceleration 
components  at  a  point  on  the  vehicle. 

3.  It  is  possible  to  represent  the  RMS  acceleration  of  a  poitit  on 
a  vehicle  as  an  integral  of  the  acceleration  power  spectrum  over  all 
frequencies.  ^  The  power  spectrum  of  a  time  dependent  or  space  depen¬ 
dent  quantity  gives  the  amount  of  power  associated  with  each  constituent 
temporal  or  spatial  frequency  component  of  the  Fourier  representation  of 
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the  physical  quantity.  ’  The  time  history  of  the  acceleration  of  a 
point  on  a  vehicle  has  a  corresponding  acceleration  power  spectrum  which 
depends  on  the  temporal  frequencies,  while  the  spatial  variations  of 
terrain  features  has  an  associated  power  spectrum  which  depends  on  the 
spatial  frequencies.  The  frequency  response  signature  for  the  accelera¬ 
tion  components  of  a  vehicle  is  the  ratio  of  the  acceleration  power  spec¬ 
trum  to  the  power  spectrum  of  the  terrain  elevation  variations.  All  of 
the  dynamical  characteristics  of  a  vehicle  are  contained  in  its  frequency 
response  signatures,  so  that  if  these  signatures  are  obtained  empirically 
they  can  serve  as  a  model  independent  representation  of  the  dynamical 
response  of  a  vehicle. 

Objectives 

4.  The  basic  objective  of  this  report  is  the  development  of  a 
semiempirical ,  model  independent,  method  of  predicting  the  dynamical 
response  of  a  vehicle  as  it  crosses  an  obstacle  or  other  terrain  feature. 
The  specific  objectives  of  this  study  are  (Figure  1): 

a.  The  employment  of  Fourier  analysis  to  describe  the  fre¬ 
quency  content  of  various  obstacle  profiles  and  the  fre¬ 
quency  content  of  the  measured  dynamic  response  of  a 
vehicle. 

_b.  The  development  of  a  method  for  determining  the  frequency 
response  signatures  of  a  vehicle  from  the  acceleration 
time  histories  measured  at  a  point  on  a  vehicle  as  it 
crosses  an  obstacle  of  known  profile. 


Scope 


5.  The  following  specific  work  was  done  (Figure  2): 

ji.  A  catalog  of  Fourier  transforms  was  developed  for  several 
basic  obstacle  shapes. 

Jb.  The  fast  Fourier  transform  (FFT)  was  used  to  determine  the 
Fourier  spectrum  of  the  acceleration  time  history  measured 
at  a  point  on  a  vehicle  as  it  crosses  an  obstacle  of  known 
profile. 

c.  Dimensionless  acceleration  frequency  response  signatures 
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Figure  2.  Scope 


were  defined  and  a  numerical  calculation  was  done  for  an 
M114  track  vehicle. 
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PART  II:  FOURIER  SPECTRUM  REPRESENTATIONS  OF  OBSTACLE 
SHAPES  AND  OF  THE  DYNAMIC  RESPONSE  OF  VEHICLES 

Introductory  Remarks 


6.  In  order  to  determine  the  frequency  response  signatures  of 
vehicles  from  obstacle  crossing  tests  it  is  necessary  to  have  both  the 
Fourier  representation  of  the  measured  dynamic  response  of  the  vehicle 
and  the  Fourier  representation  of  the  shape  of  the  cross  section  of  the 
obstacle.  The  chaotic  nature  of  the  measured  vertical  acceleration  at 
the  drivers  seat  generally  requires  that  a  computer  algorithm  be  used 
to  determine  the  Fourier  series  representation.  On  the  other  hand  the 
relatively  simple  shapes  of  obstacle  cross  sections  allows  analytical 
techniques  to  be  used  for  determining  the  Fourier  series  representation 
of  the  obstacle  geometry.  This  part  of  the  report  calculates  the  Fourier 
spectral  representations  of  the  obstacle  shapes  that  will  be  used  for 
the  determination  of  the  vehicle  frequency  response  signatures. 

7.  The  geometrical  shapes  of  actual  obstacles  expected  to  be  en¬ 

countered  on  a  battlefield  can  be  approximated  by  several  basic  geometric 
shapes — semicircle,  rectangle,  trapezoid  and  triangle.  Fourier  trans¬ 
form  and  series  representations  of  these  basic  shapes  are  easily  obtained 
by  standard  analytical  techniques.^  The  vertical  acceleration  at 

the  drivers  seat  is  generally  measured  as  a  function  of  time  as  a  vehicle 
passes  over  an  obstacle.  The  Fourier  spectrum  of  the  acceleration  time 

history  is  obtained  using  the  fast  Fourier  transform  (FFT) 

7  12  13 

algorithm.  ’  ’  The  ratio  of  the  acceleration  Fourier  spectrum  to  the 

obstacle  shape  Fourier  spectrum  essentially  determines  the  vehicle's 
frequency  response  signature. 

Fourier  Transform  Representation  of  Obstacle  Shapes 

8.  The  Fourier  transform  representation  can  be  used  to  describe 
a  single  obstacle  of  finite  extent.  Fourier  transforms  are  easier  to 
calculate  if  the  shape  of  the  obstacle  cross  section  is  symmetrical  about 


V 
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a  vertical  axis.  If  the  axis  of  evaluation  of  a  symmetrical  obstacle  is 
taken  to  be  the  axis  of  symmetry  the  Fourier  transform  reduces  to  the 
Fourier  cosine  transform  (Figure  3) .  If  the  axis  of  symmetry  of  a 


Figure  3.  Typical  Fourier  cosine  transform 


symmetrical  obstacle  is  not  chosen  as  the  axis  of  representation  or  if 
the  obstacle  is  asymmetrical,  both  Fourier  cosine  and  Fourier  sine  trans¬ 
forms  are  required.  The  Fourier  cosine  transform  pairs  are  written  as^ 


fJIoo 


cos  k£d£ 


h(£) 


cos  k£dk 


(1) 


(2) 


where 

h ( C)  *  height  of  obstacle  cross  section 

Fc(k)  ■  Fourier  cosine  transform  of  obstacle  shape 

k  =>  wave  number  *  2tt/X  where  X  =  spatial  wavelength 
components 


rife 


Jjyk 
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* 
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Defining  the  spatial  frequency  as  £2  *  1/X  yields  the  following  rela¬ 
tion  k  =  2ir£2  . 

9.  For  the  case  of  an  asymmetrical  obstacle  shape  the  following 
complex  Fourier  transform  is  used^ 


lh(k) 


hCxJe^^dx 


=  l\k)  +  ilNk) 
c  s 


Ih(k)e 


i*h(k) 


h(x)  -  —  f  Ih(k)e_ikxdk 

fa  J 


oo 

■lit 


I^(k)  cos  (kx  -  ^^)dk 


(3) 


(4) 


where 


lc<k> 


h(x)  cos  kxdx 


(5) 


I^(k) 


sin  kxdx 


Ih(k) 


M  f  +  W/ 


(6) 


(7) 
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$^(k)  =  tan 


(8) 


Ic<k). 

and  where 

I^(k)  =  complex  Fourier  transform 

l^(k)  =  real  part  of  I^(k) 
c 

I^(k)  =  imaginary  part  of  I^(k) 
s 

I^(k)  =  magnitude  of  I^(k) 

h  ~h 

$  (k)  =  phase  angle  of  the  complex  Fourier  transform  1  (k) 

10.  For  the  special  case  of  a  symmetrical  obstacle  of  length 

2L  ,  with  the  origin  of  the  axis  chosen  at  the  left  edge  of  the  obstacle 
o 

and  the  symmetry  axis  at  Lq  ,  the  Fourier  transform  components  are 
given  in  Appendix  A 


Ih(k)  =  Fh(k)  cos  kL  (9) 

c  c  o 

Ih(k)  =  Fh(k)  sin  kL  (10) 

s  c  o 


Ih(k) 


F^k) 


(ID 


4>h(k)  =  kLQ  +  nTT  (12) 

For  a  symmetrical  obstacle  with  the  origin  of  the  axis  chosen  to  be  at 
the  axis  of  symmetry  the  Fourier  transform  is  obtained  as 


ll!(k)  =  F*(k) 
c  c 


(13) 


I>)  •  0 


(1A) 
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Ih(k)  =  jF^(k) |  (15) 

$h(k)  =  nit  (16) 


which  Is  just  the  Fourier  cosine  transform. 

11.  The  dynamical  response  (acceleration)  of  several  military 
vehicles  encountering  obstacles  has  been  measured  for  basic  obstacle 
shapes  including  the  semicircle,  rectangle  and  trapezoid.  Therefore 
the  Fourier  transforms  of  several  basic  symmetrical  and  asymmetrical 
obstacles  have  been  evaluated  in  Appendix  A  and  are  tabulated  below: 

Semicircle  (Figure  4a) 


F^(k) 

1*00 


Ih(k)  =  Fh(k)  sin  kR 
s  c 


Ih(k) 


F*(k) 


Rectangle  (Figure  4b) 

-  £  4^) 


I^(k)  =  F^(k)  cos  ka 


Ih(k)  -  F^(k)  sin  ka 
s  c 


Ih(k) 


Fc'k> 
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>(R  +  s) 

2L  =  2a 

+  2s 

Semicircle 

b. 

Rectangle 

.  =  2(a  +  b  +  s) 

Symmetrical  trapezoid 


2L  =  2(t>  +  s; 

d.  Symmetrical  triangle 


c  a  a  c  c  s 

2L0  =  2(a  +  c) 

2L  =  2(a  +  c  +  s) 
e.  Double  rectangle 
Figure  4.  Symmetrical  obstacle  shapes 


Symmetrical  trapezoid  (Figure  4c) 


Fhfk>  -  i  bV.["sin  k^a  +  W2)"|/sin  kb/2\ 

Fc(k)  Vrr  \f  +  2/h[  k(r+  bjl)  JV  kb72"  / 


Ih(k)  =  Fh(k)  cos  k(a  +  b) 
c  c 


I^(k)  =  Fh(k)  sin  k(a  +  b) 
s  c 


Ih(k)  -  F*(k) 


Symmetrical  triangle  (Figure  4d) 


F^(k) 


sin  bk/2' 


i^<k) 


Fc<k)  cos  kb 


Ih(k)  =  Fh(k)  sin  kb 
s  c 


lh(k)  =  F*(k) 


Double  rectangle  (Figure  4e) 


‘(k)  -  VF  +  h2C(Sikc/f  /2)  cos  k(a  +  c/2)] 


I^(k)  ■  F^(k)  cos  k(a  +  c) 


Ig(k)  *  Fb(k)  sin  k(a  +  c) 


Ih(k)  -  Fk(k) 


Figure  5.  Asymmetrical  obstacle  shapes 
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taken  larger  than  the  width  of  the  obstacle  cross  section  (Figures  4 
and  5) .  The  Fourier  series  can  then  be  written  as^ 

A*1  °° 

h(x)  =  ~  ^  ^A^  cos  k^x  +  B^  sin  k^x\  (17) 

n=l 


where  the  Fourier  coefficients  are  given  by 


Li 

■if 


A"  =  j-  h(x)  cos  kRxdx  n  =  0,  1,  2,  3, 


(18) 


4-;| 


B“  =  r-  /  h(x)  sin  k  xdx  n=  1,  2,  3,  . . 

n 


(19) 


and  where  the  wave  number  kn  is  given  by 


nji 

L 


(20) 


The  values  of  these  coefficients  depend  on  the  position  chosen  for  the 
axis  used  to  describe  the  shape  function  h(x)  . 

13.  The  Fourier  series  can  also  be  written  in  the  following  form 


where 


h(x) 


n=l 


cos 


mix 
,  L 


(21) 


(22) 


19 


(23) 


cos 


h 


n 


h 

B  =  C  sxn  6 
n  n  n 


(24) 


(25) 


For  the  special  case  of  a  symmetrical  obstacle  cross  section  with  the 
axis  for  the  representation  of  h(x)  chosen  to  be  the  axis  of  symmetry, 
the  Fourier  cosine  series  can  be  used  as  follows"^ 


h(x) 


(26) 


where 


h 

a 


n 


mrx 
cos  — 


dx 


(27) 


for  n=0,  1,  2,  3, _  . 

14.  For  a  symmetrical  obstacle  of  length  2Lq  with  the  coordinate 
axis  chosen  at  the  left  hand  edge  of  h(x)  ,  the  Fourier  series  coeffi¬ 
cients  are  obtained  from  Equations  18,  19,  and  27  to  be 


.  h  h  .  _ 

A  =  a  cos  k  L 
n  n  no 


(28) 


B*1  =  a^  sin  k  L 
n  n  no 


(29) 


n 


n 


(30) 


20 


(31) 


<j>h  =  k  L 
n  no 

where  a||  is  given  by  Equation  27.  For  the  case  of  a  symmetrical  ob¬ 
stacle  with  the  axis  of  representation  of  h(x)  taken  to  be  the  axis  of 
symmetry,  the  Fourier  coefficients  are  given  by 


.h  h 
A  =  a 
n  n 


(32) 


Bh  =  0 
n 


(33) 


Ch  =  ah 
n  n 


(34) 


(35) 


which  is  the  Fourier  cosine  series  of  Equations  26  and  27. 

15.  A  simple  relationship  exists  between  the  Fourier  transform 
functions  and  the  Fourier  series  coefficients  for  obstacles  of  period 
2L  .  The  Fourier  series  coefficients  can  be  obtained  by  making  the  re¬ 
placements  -*■  L  and  k  -*■  k  .  With  these  replacements  it  follows 

10  n 

from  Equations  1-27  that  : 


ah  =  Fh{k  } 
n  L  c  n 


Ah  =  !h(k  } 
n  L  c  n 


Bh  .  -21  jh(k  , 

n  L  s  n 


CS  ■  T  lh(kn> 


(36) 

(37) 

(38) 

(39) 
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(40) 


h  ,  /Ih(k  ) 

ih  — 1  /  s  n 

=  tan  I  — r - 


Fourier  Spectrum  Representation  of  the 
Dynamic  Response  of  a  Vehicle 

16.  A  frequently  measured  quantity  in  vehicle  ride  quality  studies 
is  the  time  history  of  the  vertical  acceleration  produced  at  the  driver's 
seat  when  the  vehicle  encounters  an  obstacle.  This  dynamic  response  is 
generally  irregular  and  chaotic,  and  numerical  methods  must  be  used  for 

a  Fourier  analyses.  The  FFT  algorithm  is  well  suited  for  obtaining  the 
Fourier  series  coefficients  for  the  measured  acceleration  time  history. 
The  corresponding  Fourier  transforms  can  then  be  calculated  in  terms  of 
the  Fourier  series  coefficients.  The  Fourier  transforms  and  Fourier 
series  of  the  dynamic  response  will  now  be  written  out  explicitly  because 
these  quantities  appear  in  the  definition  of  the  frequency  response 
signatures  that  are  calculated  in  Part  III. 

17.  The  Fourier  transform  pairs  that  describe  the  acceleration, 
velocity,  and  displacement  at  the  driver's  seat  can  be  written  as 


a  (t) 


'a.  .  -iwt, 

I  (w)e  dw 


(41) 


v(t) 


1 

/2tT 


/ 

— CO 


~v,  .  -iwt, 

I  (w)e  dw 


(42) 


00 


;d 


(w) 


-iwt 


dw 


(43) 


where 

a(t)  =  time  history  of  the  vertical  acceleration 
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mi*  maim 


v(t)  =  time  history  of  the  velocity 
d(t)  =  time  history  of  the  displacement 
The  Fourier  transforms  are  then  written  as 


1%) 


00 

--f 

fa  J 


a(t)e*utdt 


IV(w) 


00 

fa  J 


v(t)e^utdt 


Id(ui) 


00 

■*/ 


d(t)eiutdt 


(44) 


(45) 


(46) 


where 

Ia(w)  =  complex  Fourier  transform  of  the  acceleration  at  the 
driver's  seat 

IV(u>)  =  complex  Fourier  transform  of  the  velocity  at  the  driver's 
seat 

•»  J 

1° (m)  =  complex  Fourier  transform  of  the  displacement  at  the 
driver's  seat 

oi  =  2xf  =  angular  frequency 

18.  The  complex  Fourier  transforms  can  be  written  as^ 


1%)  =  Ia(u>) 

+  ila(u>)  =  la(ui)e^a(u)) 

(47) 

IV(oi)  =  1^(01) 

v 

,  .Tv,  .  Tv,  .  i<t>  (ui) 

+  11^(10)  =  I  (w)e 

(48) 

id(o>)  =  id(to) 

.  d,  .  Td,  .  itd(ui) 

+  il  (ui)  *  I  (ui)e 
s 

(49) 
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(59 


CO 


d(t)  =  — x  +  ^  cos  a)  t  +  sin  u 

L  \  n  n  n  n  / 

n=l 

(64) 

CO 

=  2  Ct  C0S  (v  -  *n) 

n=0 


where 

u>  =  nir/T 
n 

2T  =  period  chosen  to  represent  the  time  history 
The  relationship  between  the  Fourier  series  representation  and  the 
Fourier  transform  representation  of  the  dynamic  response  of  the  vehicle 
is  given  by 


A.,V,d  .  jE  ir.<%n, 


Ba,v,d  ,  JE.  I».v.a(„n) 


,a,v,d  .  Jp  ) 


.a,v,d  _  ,a,v,d 


rn 


<T’v‘“<V 


(65) 

(66) 

(67) 

(68) 


Either  Fourier  representation  can  be  obtained  from  the  other,  so  that 
the  representation  easiest  to  obtain  from  the  measured  data  can  always 
be  used  to  determine  the  alternate  representation. 

20.  Generally  only  the  time  history  of  the  acceleration  at  the 
drivers  seat  is  measured,  but  the  velocity  and  displacement  time  histo¬ 
ries  can  be  obtained  by  numerical  integration  with  proper  choice  of 
initial  conditions.  The  measured  acceleration  time  history  is  generally 
very  complicated  and  chaotic,  and  an  accurate  determination  of  the 

26 


Fourier  transforms  by  direct  integration  is  expensive  and  time  consuming. 
Therefore  an  alternative  method  for  calculating  the  Fourier  transforms 
or  series  is  required.  A  rapid  and  inexpensive  method  of  calculating 
the  Fourier  series  coefficients  is  given  by  the  fast  Fourier  transform 
(FFT) . 


21.  The  FFT  is  a  mathematical  procedure  for  rapidly  determining 

12  13 

the  Fourier  series  of  a  specified  function.  *  The  algorithm  requires 
2n  pairs  of  data  points  t  ,  a(t  )  from  which  half  that  number  of 
Fourier  coefficient  pairs  A^  ,  are  determined.  The  accuracy  of 

the  calculated  Fourier  series  coefficients  improves  as  the  number  of  data 
points  describing  the  function  increases.  The  Fourier  transforms  can 
then  be  estimated  in  terms  of  the  Fourier  series  coefficients  by  using 
Equations  65-68. 

22.  The  vehicle  frequency  signature  that  is  calculated  in  Part  III 
is  essentially  a  comparison  (ratio)  of  the  Fourier  representation  of  a 
vehicle  response  to  the  Fourier  representation  of  an  obstacle  shape. 

These  Fourier  representations  must  be  compatible  in  the  sense  that  both 
must  be  either  Fourier  transforms  or  Fourier  series.  Also,  the  half¬ 
period  L  of  the  Fourier  series  representation  of  an  obstacle  shape 
must  be  related  to  the  half-period  T  of  the  Fourier  series  representa¬ 
tion  of  an  acceleration  time  history  by  the  equation  L  =  uT  where 

u  =  vehicle  speed  which  is  assumed  to  be  constant.  From  Equations  39 
and  67  it  follows  that 


Ca’V’d 

n 


n 


*a’v>d(y 

lh<V 


(69) 


Equation  69  is  basic  to  the  calculation  of  the  frequency  response  signa¬ 
tures  done  in  Part  III. 


PART  III:  VEHICLE  FREQUENCY  RESPONSE  SIGNATURES 


Introductory  Remarks 

23.  The  frequency  response  signature  of  a  vehicle,  which  will  be 
defined  and  calculated  here,  is  essentially  the  ratio  of  the  power  spec¬ 
trum  of  the  acceleration  response  of  a  vehicle  to  the  power  spectrum  of 
the  terrain  feature  producing  the  acceleration  response.  The  accelera¬ 
tion  response  of  a  vehicle  generally  has  three  components  that  are  asso¬ 
ciated  with  the  heaving,  pitching,  and  rolling  motion  of  the  vehicle,  so 
that  three  distinct  frequency  spectrum  signatures  are  required  for  a 
complete  description  of  each  vehicle.  Knowledge  of  the  frequency  re¬ 
sponse  signatures  makes  possible  the  prediction  of  the  vehicle  response 
to  any  terrain  features. 

24.  Model  dependent  expressions  for  the  frequency  response  signa¬ 
tures  have  already  appeared  in  the  literature,  and  in  principle  these 
signatures  are  vehicle  characteristics  that  depend  on  the  dynamical  re¬ 
sponse  parameters  of  a  vehicle  and  on  the  geometry  of  the  vehicle-ground 
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contact  area.  *  This  part  of  the  report  gives  a  model  independent 
method  of  determining  the  frequency  response  signatures  from  the  accel¬ 
eration  time  histories  measured  at  a  point  on  a  vehicle  as  it  crosses  an 
obstacle  of  known  shape. 

Calculation  of  the  Frequency  Response  Signatures 

Definitions  and  terms 

25.  The  vehicle  frequency  response  signatures  are  dimensionless 
functions  that  relate  the  dynamic  response  of  a  vehicle  to  the  tertain 
elevations  that  produce  the  response.  The  frequency  content  of  the  ter¬ 
rain  elevation  variations  can  be  described  by  a  frequency  spectrum  analy- 
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sis  using  power  spectra.  ’  ’  Obstacles  can  be  handled  in  a  similar 
fashion  and  have  the  advantage  that  their  Fourier  spectra  are  easily 
calculated.  A  simple  relationship  exists  between  the  power  spectrum  and 
the  Fourier  amplitudes  of  elevation  variations,  so  that  the  power  spectra 


associated  with  the  Fourier  spectra  of  Part  II  are  easily  calculated. 

26.  The  power  spectrum  describes  the  frequency  content  of  the 
variations  of  a  physical  quantity  in  such  a  way  that  the  integral  of  the 
power  spectrum  over  all  frequencies  is  related  to  the  root  mean  square 
(RMS)  of  the  physical  variable. ^»^>9»10  For  instance,  for  a  terrain 
elevation  such  as  an  obstacle 


1 


2 

h 


Ph(fi)dfl  = 


/ 


Ph(f)df 


(70) 


where 

E,  =  RMS  of  the  obstacle  profile 
n 

P^(fi)  =  power  spectrum  of  the  obstacle  profile  expressed  in  terms 
of  the  spatial  frequency 

P^(f)  -  power  spectrum  of  the  obstacle  profile  expressed  in  terms 
of  the  time  frequency 

and  where  the  time  frequency  is  related  to  the  spatial  frequency  by 

f  =  uft  (71) 

where  u  equals  the  horizontal  speed  of  the  vehicle  which  is  assumed  to 
be  constant  during  transit  of  the  obstacle.  From  Equations  70  and  71 
it  follows  that 

Ph(n)  =  uPh(f)  (72) 

27.  A  complete  description  of  the  dynamic  response  of  a  vehicle 
includes  the  time  histories  of  the  acceleration,  velocity,  and  displace¬ 
ment  of  a  point  on  the  vehicle  during  transit  of  the  obstacle.  In  prac¬ 
tice  only  the  three  components  of  the  acceleration  are  measured.  How¬ 
ever,  for  completeness  the  frequency  response  signatures  of  the 
acceleration,  velocity,  and  displacement  will  be  considered.  The  fre¬ 
quency  spectrum  signatures  relate  the  power  spectra  of  the  dynamic 
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response  of  the  vehicle  to  the  power  spectrum  of  the  obstacle  profile 
and  are  defined  as  follows: 


CO  00 

l2a  =  /pa(f)df  =  J  (27if)4Sa(f)Ph(f)df 


(73) 


00  00 

zl  =  J  Pv(f)df  =  J  (21Tf)2Sv(f)Ph(f)df 


(74) 


00  oo 

t;  ■/ pd(f)d' s  /  sd(t>ph(f>d£ 


(75) 


where 


l 


v 


P  (f) 

%«n 

pd(f) 
S  (f) 

Vo 

sd<0 


RMS  of  the  acceleration  of  a  point  on  the  vehicle  for  one 
ot  the  three  inodes  of  motion 

RMS  of  the  velocity  of  a  point  on  the  vehicle  for  one  of 
the  three  possible  modes  of  motion 

RMS  of  the  displacement  of  a  point  on  the  vehicle  for  one 
of  the  three  possible  modes  of  motion 

power  spectrum  of  the  acceleration  modes 

power  spectrum  of  the  velocity  modes 

power  spectrum  of  the  displacement  modes 

acceleration  frequency  response  signatures 

velocity  frequency  response  signatures 

displacement  frequency  response  signatures 


Because  the  acceleration,  velocity,  and  displacement  each  has  three  com¬ 
ponents,  there  are  in  principle  nine  frequency  response  signatures.  In 
practice  only  the  three  acceleration  frequency  response  signatures  are 
measured. 


Frequency  response  signa¬ 
tures  for  acceleration, 
velocity  and  displacement 

28.  From  Equations  72-75  it  follows  that  the  frequency  response 
signatures  are  given  by 

P .(f)  uP  (f) 

s  (f)  =  - a_ - * -  (76 

(2irf )  Ph(f)  (2nf)UP  (fl) 


P  (f)  UP  (f) 

Vf) - ~2 -  =  - 2 - 

(2rf)ZP  (f)  (2*0^  (fi) 


P  (f)  uP  (f) 

Q  (f)  -  - -  =  - 9 - 

bd(t)  Ph(f)  Ph(C2) 


For  the  purpose  of  numerical  calculations  it  is  necessary  to  express  the 
power  spectra  in  Equations  76-78  in  terms  of  Fourier  amplitudes  which 
are  obtained  from  a  frequency  spectrum  analysis  of  the  obstacle  profile 
and  the  measured  time  histories  of  the  dynamical  vehicle  response. 

29.  The  power  spectrum  of  a  function  is  related  to  the  square  of 
the  magnitude  of  the  Fourier  transform  of  the  function  as  follows 


p  (f)  =  f  !ia(u») 
a  o 


P  (f)  =  f  |lV(u>) 


P,(f)  =  f  |l%) 
d  o 


Ph(n)  =  (k) 


where 


u>  *  2irf 
k  -  2rJ2 


The  Fourier  transforms  were  calculated  in  Part  II.  The  constants  f 

o 

and  are  chosen  for  normalization  purposes  and  are  related  by 

f  =  ufl  corresponding  to  Equation  71.  The  value  of  f  is  taken  from 
O  u  o 

the  measured  duration  of  the  time  history  2T  as  f^  =  1/(2T)  . 

30.  From  Equations  76-82  it  follows  that  the  frequency  response 
signatures  can  be  written  as 


S  (u>) 
a 


Sv(u>) 


Sd(-) 


(83) 


(84) 


(85) 


where  u>  and  k  are  releated  by  <*>  *  uk  .  Therefore  the  dimensionless 
frequency  response  signatures  can  be  expressed  in  terms  of  the  ratios 
of  the  Fourier  transforms  of  the  vehicle  response  to  the  Fourier  trans¬ 
form  of  the  obstacle  profile. 

31.  The  frequency  response  signatures  can  be  rewritten  in  terms 
of  transmission  functions  which  are  defined  as  follows 


S  (to)  =  tJ(u) 
a  a 

(86) 

/■s 

3 

v 

JJ 

3 

> 

CO 

(87) 

Sd(w)  =  tJ(io) 

(88) 

where 

T  (w)  =  acceleration  transmission  function 
a 

Tv(w)  =  velocity  transmission  function 


32 


T^Cto)  *  displacement  transmission  function 
which  are  found  from  Equations  83-85  to  be 


t  (»>  -  4 

u2  Ih(k) 


I  (ii))  *  —  V551 
v  “  Ih(k) 


T,(i>)  ■  "V1^  (91) 

d  l\k) 

32.  Because  the  dynamic  response  time  histories  of  a  point  on  a 
vehicle  are  very  complex  it  is  more  convenient  to  use  the  Fourier  series 
coefficients  as  calculated  in  Part  II  by  the  FFT  rather  than  the  Fourier 
transforms  to  calculate  the  frequency  response  signatures.  The  connec¬ 
tion  between  the  Fourier  series  coefficients  and  the  Fourier  transforms 
is  given  by  Equation  69  and  when  this  is  combined  with  Equations  83-91 

and  a)  =  mr/T  the  following  results  are  obtained 
n 


■  d?)*  g 


V”>  ■  (hf  (| 


=  T2(n) 
a 


=  T2  (n) 
v 


Sd(n)  = 


where 


T  (n) 
a 


-jr  ’  Ta(n) 

a/  d 


(T-)2C- 
VW  „h 


Tv(n) 


IMI 

n 


V"> 


for  n=l,  2,  3,  ...  .  For  n  =  o  the  frequency  response  signatures 
are  by  definition  equal  to 


Sfl(o)  =  Sv(o)  =  Sd(o)  =  1 


Ia(0)  -  Tv(o)  -  Td(o)  -  1 


Table  1  gives  the  dimensions  of  all  the  physical  quantities  that  appear 
in  this  report. 

33.  From  Equations  62-64  and  w  =  rnr/T  it  follows  that 

n 


T  (n)  =  T  (n) 
a  v 


Td<„> 


WOO) 


This  follows  from  the  assumption  that  the  acceleration,  velocity  ahd 
displacement  of  a  point  on  the  vehicle  can  all  be  expanded  as  Fourier 
series  with  a  common  period  2T  as  shown  in  Equations  62-64. 

34.  Having  obtained  Sa  ,  and  from  response  data  for  a 

vehicle  on  a  test  obstacle,  it  is  now  possible  to  predict  the  RMS  values 
of  acceleration,  velocity  and  displacement  that  are  expected  to  occur 
at  a  point  in  the  vehicle  when  it  crosses  another  obstacle  of  known  pro¬ 
file.  This  is  done  using  Equations  73-75  where  now  the  power  spectrum 
P^(f)  *  P^CSty/u  refers  to  the  new  obstacle.  The  new  power  spectrum  is 
obtained  using  Equation  82  and  the  Fourier  transform  of  the  new  obsta¬ 
cle  profile.  In  this  way  the  RMS  dynamical  response  of  a  vehicle  can  be 


predicted  for  any  obstacle  or  any  terrain  elevation  variation  without 
having  a  specific  dynamical  model  for  the  vehicle.  Only  the  empirical 
frequency  response  signatures  of  a  vehicle  are  required  to  predict  the 
RMS  dynamical  response  on  any  type  of  terrain.  This  is  true  only  within 
the  limits  of  the  validity  of  the  assumption  that  the  acceleration, 
velocity,  and  displacement  of  a  point  on  the  vehicle  can  be  written  as 
Fourier  series  in  the  form  given  in  Equations  62-64. 

Model  Description  for  the  Frequency 
Response  Signatures 

35.  In  order  to  obtain  a  better  understanding  of  the  meaning  of 
the  frequency  response  signatures,  their  values  for  a  simple  vehicle- 
ground  contact  model  is  presented.  Theoretical  model-dependent  expres¬ 
sions  for  the  frequency  spectrum  signatures  have  already  appeared  in  the 
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literature  and  are  as  follows  ’  : 


-*( 


S  •  T.J  e 
a  1  dd 1 


S  =  | T  I 2  e~^  ^ 
v  1  Xdd 1  8 


s  =  It  | ^  e~F 
bd  1  dd1 


asw<4> 

wheels\ 

asw(4) 

track  / 

asw<2> 

wheels\ 

(2) 

asw 

track  / 

wheels 

or  track 

(101) 


(102) 


(103) 


where 


|  T^^ I  =  magnitude  of  the  displacement-displacement  transmission 
function  for  vehicle 


-F 


(A) 


=  low  pass  filter  associated  with  the  contact  length  of  a 
wheel  or  track 

*  curvature  spectral  window  function  for  wheeled  vehicle 


ASW 

asw^^  =  curvature  spectral  window  function  for  track 


(2) 

ASW  =  slope  spectral  window  function  for  wheeled  vehicle 

(2) 

asw  =  slope  spectral  window  function  for  track 

36.  Equations  101-103  show  that  in  principle  the  frequency  re¬ 
sponse  signatures  are  independent  of  the  obstacle  profile  and  depend 
only  on  the  geometric  and  dynamic  characteristics  of  the  vehicle.  The 
dynamic  characteristics  enter  through  the  displacement-displacement 
transmission  function,  while  the  geometry  of  the  vehicle-obstacle  con¬ 
tact  area  enters  through  the  low  pass  filter  associated  with  the  spectral 

window  functions  and  the  low  pass  filter  associated  with  the  ground- 
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vehicle  contact  length.  ’  For  a  very  small  vehicle-ground  contact 

length  the  frequency  response  spectra  in  Equations  101-103  reduce  to 
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the  square  of  the  dynamic  transmission  function.  ’ 

37.  The  contact  length  between  the  vehicle  and  the  obstacle  is 
very  small  (about  12  in.*)  for  the  case  of  a  vehicle  crossing  the  small¬ 
sized  obstacles  that  were  used  to  obtain  the  experimental  acceleration 
time  histories  treated  in  this  report.  This  is  true  for  both  track  and 
wheeled  vehicles  because  the  obstacles  used  for  the  tests  are  about 

12  in.  wide.  For  this  case  it  follows  from  Equations  101-103  that 


S  =  S  =  S  =  T.J 
a  v  d  '  dd 1 


(104) 


and  the  conditions  for  the  validity  of  Equations  100  and  62-64  are 
satisfied.  It  follows  that 

Ta  '  lTddl  (105) 

and  the  empirical  values  of  the  effective  transmission  function  should 
be  an  accurate  measure  of  the  displacement-displacement  transmission 


function  for  a  vehicle. 


Numerical  Results 


38.  A  numerical  example  of  the  theoretical  procedures  developed 
in  this  report  will  now  be  presented.  The  vertical  acceleration  time 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measurement 
to  metric  (SI)  units  is  presented  on  page  3. 
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history  at  the  drivers  seat  of  an  M114  track  vehicle  crossing  a  6-in. 

.  radius  semicircular  obstacle  at  a  speed  u  =  5.2  mph  =  91.5  in. /sec 

* 

will  be  used  to  determine  the  transmission  function  from  which  the  fre¬ 
quency  response  signature  follows  trivially  as  its  square.  The  time 
history  of  the  vertical  acceleration  at  the  drivers  seat  appears  in 
Figure  6. 

39.  The  Fourier  series  representation  of  the  acceleration  time 
history  was  obtained  from  a  measured  period  of  2T  =  1.55  sec  using 

an  FFT  computer  program  and  appears  in  Figure  7.  The  Fourier  series 
representation  of  a  6- in.  radius  semicircular  obstacle  was  obtained 

for  a  period  2L  =  2yT  =  141.825  in.  using  the  analytical  expression 
for  the  Fourier  transform  that  appears  in  the  table  of  transforms  com¬ 
bined  with  Equation  39,  and  appears  in  Figure  8.  The  transmission  func¬ 
tion  for  the  vehicle  was  then  calculated  using  Equation  95  and  the 
results  appear  in  Figure  9. 
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semicircle  obstacle 


Fourier  series  coefficients  for  acceleration 


PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

40.  This  report  develops  a  formalism  for  calculating  the  fre¬ 
quency  response  signatures  of  a  vehicle  from  the  acceleration  time 
history  that  is  measured  at  a  point  on  a  vehicle  as  it  crosses  an  obsta¬ 
cle  of  known  profile.  The  formalism  uses  the  Fourier  series  representa¬ 
tions  of  both  the  obstacle  shape  and  the  measured  acceleration  time 
history. 

41.  The  study  of  the  dynamic  response  of  vehicles  to  obstacles 
produced  the  following  conclusions: 

a.  A  nondimensional  vehicle  characteristic  called  the  fre¬ 
quency  response  signature  can  be  defined  and  calculated 
as  being  proportional  to  the  square  of  the  ratio  of  the 
Fourier  spectrum  of  the  measured  acceleration  to  the 
Fourier  spectrum  of  the  shape  of  the  obstacle. 

b.  The  acceleration  frequency  response  signatures  are  a 
measure  of  three  distinct  effects  (1)  the  dynamic  trans¬ 
mission  function  which  includes  the  effects  of  the  mass, 
elastic  spring  and  damping  parameters  of  the  vehicle 

(2)  the  low  pass  filtering  effect  of  the  track  or  wheel 
ground-contact  length  (3)  the  low  pass  filtering  effect 
of  the  curvature  spectral  window  function  for  the  vehicle 
which  accounts  for  the  curvature  of  the  obstacle  over 
the  contact  length. 

£.  A  brief  numerical  study  showed  that  the  characteristic 
resonance  peak  of  the  vehicle  displacement-displacement 
transmission  function  can  be  determined  from  an  accelera¬ 
tion  time  history  measured  at  a  point  on  a  vehicle  as  it 
crosses  an  obstacle  of  known  shape. 

Recommendations 


42.  The  study  of  the  empirical  determination  of  the  frequency 
response  signatures  for  vehicles  produced  several  recommendations  that 
will  be  helpful  for  the  design  of  new  vehicles.  It  is  recommended 
that  the  following  work  be  done: 

a.  Frequency  response  signatures  should  be  determined 


for  a  series  of  wheel  and  tracked  vehicles  and  used 
to  extract  the  effective  damping  and  spring  constants 
for  these  vehicles. 

b.  The  frequency  response  signatures  should  be  deter¬ 
mined  for  a  series  of  vehicle  speeds  in  order  to 
examine  the  effects  of  nonlinear  suspension  systems. 
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Table  1 

Dimensions  of  Physical  Quantities 


Units 

Obstacle  profile  height  h(x) 

In. 

Displacement  at  a  point  on  vehicle  d(t) 

in. 

Velocity  at  a  point  on  vehicle  v(t) 

in. /sec 

Acceleration  at  a  point  on  vehicle  a(t) 

in . / sec 

Fourier  transform  of  obstacle  profile  I*1 

00 

in.2 

Fourier  transform  of  displacement  I^(f) 

in .  sec 

Fourier  transform  of  velocity  IV(f) 

in. 

Fourier  transform  of  acceleration  Ia(f) 

in. /sec 

Fourier  coefficients  of  obstacle  profile 

ch 

n 

in. 

Fourier  coefficients  of  displacement 

in. 

y 

Fourier  coefficients  of  velocity  Cn 

in . / sec 

Fourier  coefficients  of  acceleration  C 

XI 

in. /sec' 

Power  spectrum  of  obstacle  profile  P^Cfl) 

in.3 

Power  spectrum  of  displacement  P^(f) 

2 

in .  sec 

Power  spectrum  of  velocity  Pv(0 

in.2/sec 

Power  spectrum  of  acceleration  P  (f) 

in.2/sec 

Time  period  of  acceleration  2T 

sec 

Space  period  of  obstacle  2L  *  2uT 

in. 

Obstacle  length  2L 

in. 

APPENDIX  A:  CALCULATION  OF  FOURIER  TRANSFORMS 
OF  OBSTACLE  SHAPES 


1.  This  appendix  gives  a  sketch  of  the  calculations  used  to  deter¬ 
mine  the  Fourier  transforms  of  the  basic  obstacle  shapes  considered  in 
this  report.  The  problem  is  to  calculate  the  Fourier  cosine  and  sine 
transforms  defined  as 


c 


cos  kxdx 


(Al) 


I 


h 

s 


sin  kxdx 


(A 2) 


where  h(x)  equals  the  height  profile  of  the  obstacle.  The  vehicle  is 
assumed  to  be  moving  from  left  to  right,  so  that  the  left  edge  of  the 
obstacle  is  encountered  first  and  it  is  at  this  point  that  the  accelera¬ 
tion  response  time  history  of  the  vehicle  begins.  Therefore  the  origin 
of  the  axis  is  placed  at  the  left  edge  of  the  obstacle.  For  an  obstacle 
of  length  2Lq  the  integrals  in  Equations  Al  and  A2  can  be  written  as 

2L 

o 

I*1  =  _1_  /  h(x)  cos  kxdx  (A3) 

c  »/2tT  J 


I 


h 

s 


h(x)  sin  kxdx 


(A4) 


where  x  equals  the  coordinate  measured  from  the  left  edge  of  the 
obstacle. 


2.  It  Is  helpful  to  introduce  the  following  coordinate 
transformation 


x  =  L  +  C 
o 


Then  the  integrals  in  Equations  A3  and  A4  can  be  written  as 


s 


h(C)(cos  kLQ  cos  k£  -  sin  kLQ  sin  k£)d£ 


h(£)(sin  kLQ  cos  k?  +  cos  kLo  sin  k£)d£ 


If  the  obstacle  is  symmetrical  the  sine  integrals  in  Equations  A6 
vanish  and  the  cosine  integrals  are  symmetric  so  tha-.. 


c 


cos  k£d£ 


cos  k£d£ 


which  can  be  rewritten  as 


c 


c 


cos  kL 

o 


s 


sin  kL 
C  0 


(A5) 

(A6) 

(A7) 

and  A7 

(A8) 

(A9) 

(A10) 

(All) 


A2 


where  F  is  the  Fourier  cosine  transform  defined  as 
c 


c 


h(£)  cos  k£d£ 


(A12) 


Therefore  for  symmetrical  obstacles  it  is  only  necessary  to  calculate 

tl  ll  tl 

the  Fourier  cosine  tranform  F  because  I  and  I  follow  automati- 

c  c  s  , 

cally  from  Equations  A10  and  All.  For  asymmetrical  obstacles  both  i” 
and  1^  must  be  calculated  separately. 

3.  The  following  parts  of  this  appendix  gives  the  evaluation  of 
the  Fourier  transforms  of  some  basic  obstacle  shapes.  Symmetrical  obsta¬ 
cles  are  treated  first  and  then  their  asymmetrical  counterparts  are 
considered . 


Symmetrical  Obstacles 


Semicircular  obstacle  (Figure  4a) 

4.  The  semicircle  is  described  by  Lq  «*  R  and 

h(£)  =  R>/l  -  £2/R2  (A13) 

for  |S  |  £R  ,  where  S  is  measured  from  the  symmetry  axis.  The  func¬ 
tion  is  symmetrical  and  the  Fourier  cosine  transform  is  given  as 


F^(k) 

c 


cos  k£dS 


(A14) 


The  substitution  n  -  S/R  gives 


F*<k) 


cos  (kRn)dn 


(A15) 


But  from  the  table  in  Sneddon's  book^  it  follows  that 


VI/ 


4/  VTT7 


cos  (kRn)dn 


J, (kR) 


nr  il 

\  2  kR 


(A16) 


where  J^(kR)  equals  the  first  order  Bessel  function.  Then  it  follows 
that  the  Fourier  transform  of  the  semicircle  is 

J, (kR) 


Fh(k)  =  R2 
c 


VP 


(A17) 


For  small  kR  one  has  J^(kR)  <v.  kR/2  ,  so  that 


Ak) 

c 


1  [F 

2  \  2 


R2  =  0.62666R2 


(A18) 


For  large  kR 


Jj.CkR)  ^  C0S  (kR  "  3lT/2) 


(A19) 


and  the  Fourier  transform  is  obtained  from  Equation  A17  to  be 


■>>  *  Vf  r  VS cos  <kR ' 3,/2) 


(A20) 


The  zeros  of  J^kR)  are  given  in  the  references  and  are  determined  by 


Jl(knR)  =  0 


(A21) 


which  gives 


k  R 
n 


r 

n 


A4 


(A22) 


where  r  equals  the  n'th  root  of  the  first  order  Bessel  function.  The 
n 

roots  of  are  tabulated  in  the  literature  and  the  first  few  are  as 

follows 


Therefore  the  values  of  k 


3.8317 

(A23) 

7.0156 

(A24) 

10.1735 

(A25) 

and  X  for  the  roots  are 

r 

II 

R 

(A26) 

r 

_  II 

2xR 

(A27) 

=  2xR 

V 

(A28) 

M, 

n 

Rectangular  obstacle  (Figure  4b) 

5.  The  rectangle  is  described  by  h(g)  =  h  for  0  <  5  <  a  where 
£  is  measured  from  the  axis  as  symmetry,  and  the  Fourier  cosine  trans¬ 
form  for  this  symmetrical  obstacle  with  Lq  =  a  is 


(A29) 


The  zeros  occur  at 


(A30) 


k  a  «  tin 
n 


The  values  of  k  , 


ft  ,  and  X  at  the  roots  are 


.  nit 

kn  =  “a 

(A31) 

or 

ft  *  f- 
n  2a 

(A32) 

X  =2a 
n  n 

(A33) 

Symmetrical  trapezoid  (Figure  4c) 

6.  The  symmetrical  trapezoid  is  described  by  the  following  func¬ 
tions  with  the  origin  of  coordinates  chosen  at  the  axis  of  symmetry  and 

L  =  a  +  b 
o 

(h  0  <  c  <  a 

(A34) 

£  (a  +  b  -  K)  a  <  £  <  a  +  b 

where  the  coordinate  £  is  measured  from  the  symmetry  axis.  The 
Fourier  cosine  transform  is  given  by 


F^k) 


cos  k£d£  + 


£  (a  +  b  -  O  cos  k?d£ 


=  ^2  [cos 

bk 


ka  -  cos  k(a  +  b)] 

tsin  k(a  +  b/2)~|  [cin  (kb/ 2)1 

k(a  +  b/2)  J  L  kb/2  J 


(A35) 


The  limit  b  -*■  0  reproduces  the  results  for  the  rectangle  shown  in 
case  2. 


A6 


& 


Symmetrical  triangle  (Figure  4d) 

7.  The  case  of  the  symmetrical  triangle  can  be  obtained  from  the 
symmetrical  trapezoid  by  taking  a  ■  0  .  From  the  result  in  case  3  for 
the  symmetrical  trapezoid  it  follows  that  the  symmetrical  triangle  is 
obtained  with  a  ■  0  to  be 

and  the  obstacle  half-length  is  Lq  =  b  . 

Double  rectangle  (Figure  4e) 

8.  The  double  rectangle  obstacle  can  be  represented  as  follows 


h(C)  - 


0  <  £  <  a 
a  <  £  <  a  +  c 


(A37) 


where 

£  ■  coordinate  measured  from  the  symmetry  axis 

L  -  a  +  c 
o 

The  Fourier  cosine  transform  is 


F^k) 


vi j  \  *°° >*« + v?  7 

*a 


cos  k£d$ 


-  VF  "i 

, 


sin  ka  .  ,  I  sin  k(a  +  c)  sin  ka 

__  +  h2  ^  -  j— 


i! 


(A38) 


>/![v 


sin  ka  .  v  sin  kc/2  .  _  ,„\ 

— kl~  +  h2c  '  W2“-  COS  k(a  +  c/2)J 


For  c  -  0  this  reduces  to  case  2  of  the  rectangle  of  height  h^  , 
while  for  a  -  0  it  reduces  to  case  2  of  the  rectangle  of  height  hj  . 


A7 


Asymmetrical  Obstaclea 


Asymmetrical  trapezoid  (Figure  5a) 

9.  The  origin  of  the  axis  is  taken  at  the  left  edge  of  the  obsta¬ 
cle  whose  total  length  is  2Lq  ■  2a  +  b  +  e  .  The  equations  defining 
the  shape  of  the  obstacle  are 


a  <  x  <  b 

b  <  x  <  2a  +  b  (A39) 

2a  +  b<x<2a  +  b  +  e 


The  cosine  component  of  the  Fourier  transform  is 


b  2a+b 

I*1  «  — —  /  ^  x  cos  kx  +  — /  h  cos  kxdx 

c  fa  J  b  SU  jj, 


2a+b+e 

-i-  f  i 
X  e 


(2a  +  b  +  e  -  x)  cos  kxdx 


(A40) 


After  some  algebra  the  cosine  component  is  found  to  be 

_h  1  h  F  ,  ...  e\^sin  ke/2\  kb  ^sin  kb/2^  | 

ic  ■  —  t  k(2a  +  b  +  ijl'-k?/— )  - 1 V  wj  )  J  (M1) 

The  sine  component  of  the  Fourier  transform  is  given  by 


.A 

/ui  b 


2a+b 


sin  kxdx  +  — —  /  h  sin  kxdx 

i/2v  J. 

O  D 


2a+b+e 

4-  —  /  -  (2a  +  b  +  e  - 

*  zL  e 


(A42) 


x)  sin  kxdx 


A8 


Some  algebra  gives  the  sine  component  of  the  Fourier  transform  as 


^rk  L 


sin  kb  / sin  ke/2 \ 

kb  \  ke/2  ) 


cos  k^2a  +  b  + 


(A43) 


The  magnitude  of  the  Fourier  transform  is  then  obtained  as 


lh  ■  >/('«")’  ♦  (■!)' 

-  ±  i  J(w)2  *  w  -  KtetP)  W)  ~ 


(A44) 


Note  Equation  A44  is  symmetric  under  interchange  of  the  variables 
b  and  e  . 

Right  angle  trapezoid  (Figure  5b) 

10.  The  right  angle  trapezoid  can  be  obtained  as  a  special  case 
of  the  general  asymmetrical  trapezoid  by  taking  e  ■  0  in  case  1. 

The  results  are 


1  h  T  v  kb  /sin  kb/2\2 

—  ^  |sin  k(2a  +  b)  -  -j  \  kb?2  /  J 


1  h  [: 

JiH  k  L 


sin  kb 
kb 


-  cos  k(2a  +  b) 


] 


(A45) 


(A46) 


1 


h  L  .  /sin  kb/2\2 
k  \  1  +  V  kb/2  ; 


-  Ktett1)  cos  k(2‘ + !)  (M7) 


Asymmetrical  triangle  (Figure  5c) 

11.  The  asymmetrical  triangle  can  be  obtained  as  a  special  case 
of  the  general  asymmetrical  trapezoid  by  taking  a  ■  0  in  the  results 
of  case  1.  The  results  are 


Right  triangle  (Figure  5d) 

12.  The  Fourier  transform  of  the  right  triangle  can  be  obtained 
as  a  special  case  of  the  asymmetrical  triangle  (case  3)  by  taking  either 
b  =  0  or  e  =  0  depending  on  whether  the  right  angle  faces  the  moving 
vehicle  or  not. 

a..  Right  angle  facing  vehicle  (b  =  0) 

.2 


h  h  /ke\/sin  ke/2\ 


ih- 

s 


/2iTk 


sin  ke\ 
ke  / 


I"  = 


/2irk 


sin  ke/2\ 
ke/2  / 


-  2 


sin  ke 
ke 


b.  Right  angle  facing  away  from  vehicle  (e  *=  0) 

kb  /sin  kb/2\2l 


(A31) 


(A52) 


(A53) 
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